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Abstract 

We introduce a method of constructing a general Laakso space while 
calculating the spectrum and multiplicities of the Laplacian operator on 
it. Using this information, we find the leading term of the trace of the 
heat kernel and the spectral dimension on an arbitrary Laakso space. 
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1 Introduction 



Much work has been done on the analysis of fractals, specifically concentrating 
on the spectrum of the Laplacian operator on irregular domains. One such 
topic are drums with Koch snowflake boundary, see for example 13J. This 
paper will be concerned instead with the irregular domain being a fractal itself. 
Some notable works with this type of domain include [31 [H [HI HE1 El US] among 
others. Laakso's spaces were introduced in [TTJ. They are a family of fractals 
with an arbitrary Hausdorff dimension greater than one and were considered 
originally for their nice analytic properties. Constructions of the Laakso spaces 
are given in [TTJ [THl [T7] as well as in Section [2] of this paper. Theorem 6.1 in 
|16j gives the spectrum of the Laplacian operator on any given Laakso space, in 
Theorem |3.1| we give the multiplicities. 

An important part of the analysis of Laplacians is the heat equation and 
associated heat kernel, which can reveal significant information about the oper- 
ator and underlying space. The information gained from studying heat kernels 
can be applied in other areas of analysis as well as other fields such as physics. 
The papers [5[ Qj|] are devoted to finding and analyzing the heat kernel and the 
trace of the heat kernel. The notion of complex valued fractal dimensions and 
the accompanying oscillating behavior of the heat kernel were studied in [3] . 

We begin by reviewing the construction of the Laakso spaces as presented 
in [TTJ [2)1 [T7] in Section [2j This section also contains background information 
on the Hausdorff dimension, its calculation for Laakso spaces, and some specific 
values for certain Laakso spaces. In subsection 2.3 we define the Laplacian 
operator that will be used throughout the rest of this paper. 

In Section [3] we begin by stating the spectrum of the Laplacian and the 
multiplicities of each eigenvalue (Theorem 3.1 1, while the rest of the section is 
devoted to the proof of this result. In Section |37T] we provide an analytical proof 
by examining, as in [16) . the different "shapes" that make up the space. Since 
each shape has a unique contribution to the spectrum counting the number of 
shapes allows us to calculate the spectrum with multiplicities. In Section |3.2 
we verify the results computationally using MATLAB. Finally in Sections [4 
and [5] we use the spectrum and multiplicities obtained in Sections |2.2| and |3.1 
to calculate the trace of the heat kernel using the same method outlined for 
diamond fractals in [2]. 



2 Laakso Spaces 

The spaces that will be analyzed were first defined by Laakso in [TT]. Laakso's 
spaces form an uncountable family of metric-measure spaces indexed by se- 
quences {j'n}nLi- An equivalent construction using projective limits was hinted 
at in [5] and fully developed in [17j . Then in 16 , a more in depth description 
of the projective limit construction was used to calculate the spectrum of the 
Laplacian constructed in [T7]. We will be using the construction in [TB] as it is 
also well-suited for our calculations. 
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The Laakso space can be visualized with a sequence of quantum graphs, 
denoted F n ,n > 0, each an increasingly better approximation of the Laakso 
space. The first of these graphs is the unit interval, denoted F . Laakso spaces 
are defined by a sequence {j n }^Lx of integers j n > 2, where each j n described 
the number of identifications at step n of the construction. To construct the 
graph of F n+ i, first every cell, or interval between two nodes, of the F n graph 
is split evenly into j n segments by adding nodes. This graph is then duplicated 
and connected at the newly-added notes. In this visualization, all nodes are 
arranged in columns. 





Figure 1: Construction of Fx from Fq with ji = 2. 

We describe a simple case, where j n = 2, n > 1. To obtain Fx, bisect F 
with a node. Then make a copy of this graph. Identify the new nodes "glueing" 
the two graphs together, represented by the arrow in Figure [TJ This glueing 
process is the identification process described in [5]. This yields the graph Fx, 
an X-shape with five nodes as seen in Figure [I] 

This procedure is repeated to obtain F 2 from F±. Nodes bisect each cell of 
Fx as seen in Figure [2] A duplicate copy of Fx is created and the two graphs are 
"glued" together at the newly added nodes. This is shown in Figure [2] where 
the solid line represents Fx and the dashed lines represent the copy of Fx. The 
j„ = 2 Laakso space is the projective limit of the sequence of graphs {Fnj^-Q 
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all produced in this manner. 




Figure 2: Construction of Fi from F\ where 22 = 2. The dashed lines represent 
the second copy of F\ with the added nodes. 

As another simple example consider j n = 3, for all n € N. Again starting 
with the unit interval, Fq, F± is constructed by splitting Fq into three subinter- 
vals and placing a node between each interval as shown in Figure [3] This graph 
is duplicated and is glued to the original graph at the newly added nodes. The 
two nodes in the middle of the figure are connected by the middle interval and 
its copy, thus creating a loop shape that is not seen in the case where j = 2. The 
outer thirds of the figure create a "V" shape, also seen in the j = 2 construction. 
These shapes, loop and "V", will be two of those considered in Section [3j 




F F 1 

Figure 3: Construction of F\ from Fq with j\ = 3. 

In this paper, we deal with the general case where j n may vary at each ap- 
proximation level n. The sequence {j n }^Li may be a constant integer, as seen 
in the previous examples. Or the sequence may alternate regularly between 
two integers. Figure [I] shows the construction of F 2 and F 3 when {j n }^i = 
{2,3,2,3,...}. The sequence {j n }%Li could even be a completely random se- 
quence of integers. In any case, it is {j n }%Li which defines the Laakso space 
and from which the properties are derived. 
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J2=[2,3] 



is = [2,3,2] 



Figure 4: Constructions of the Laakso space for {j n }^Li = {2, 3, 2, 3, ...} 
2.1 Cell Structure of a Laakso Space 

Recall that as an inverse limit system the pair [L, {F n }) come with continuous 
projection : L — > F n . 

Definition 2.1. The cell structure in a Laakso space, L, is determined by the 
pre-images under the map $„ of the cells in the graph F n which approximates 
L 

Given a space as defined by Laakso in [TT] , and the construction of the space 
as given in |16j . and the level of approximation, n, the cell structure has specific 
properties, including number of cells, N n , and the interval length, I" 1 . Both 
the number of cells and the interval length are dependent on the choice of ji for 
all i < n. 

Proposition 2.1. Each cell in F n has metric diameter 

n 

In = Y[ji (2.1) 

i=l 

where Jq = 1. In addition the number of cells is 

n 

N n = 2 n l[. h . (2.2) 

i=l 

Proof. In Fo there is a single cell, the unit interval, with metric diameter equal 
to 1. At each step in the construction the diameter of each cell in F n is j^ 1 
times that of a cell in F n -\. By induction the diameter of the cells in F n is J" 1 . 

There is a single cell in Fo. At each step of the construction there are 2 x j n 
cells in F n for every cell in F n _i. Thus the number of cells in F n is 2 n I n . □ 
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2.2 Hausdorff dimension of the Laakso Space 



In order to discuss the Hausdorff dimension of Laakso spaces we fix our choice 
of metric and measure. We use the path length metric. The measure used 
is the probability measure that gives equal mass to all cells of a given depth. 
Implicitly in the given construction, we have restricted the Hausdorff dimension 
to 1 < Q < 2. In [TB] the Hausdorff dimension Q of the Laakso Space associated 
with a constant j n at every level n is shown to be 

Here we give the Hausdorff dimension of a Laakso space associated with a general 
sequence {j n }%Li- The measure used in calculating the Hausdorff dimension is 
the projective limit of Lebesgue measure on F n scaled to have total mass one 
for all n. 

Lemma 2.1. Given sequence the Hausdorff Dimension, Q, of the cor- 

responding Laakso space is given by 



log 2" ]J 



Ji 



h = lim x , - 1 r ' = lim " / = lim 1 + a ) ' (2.3) 

n->oo I n \ n->oo log(l n ) n->oo log(l n ) 



log(2 n I n ) y 1 , log(2 n ) 




if the limit exists. 

Proof. Laakso spaces are lacunary self-similar sets as defined in [7] where the 
contraction ratios at any n are equal. The number of identifications for each 
cell at the z'th iteration is j% an d the formula that Igudesman gives in [7] can 
be given in terms of n and ji. This formula uses the number of cells, N n and 
the cell diameter, which is simply I^ 1 . In the geodesic metric, the cell length is 
also the diameter of the cell. The resulting formula is given above. □ 

While there are many sequences {j n } for which Q will not exist it is more 
relevant to our interests that for every Q there exist sequence {j n } yielding 
a Hausdorff dimension of Q. Different sequences {j n } can yield the same di- 
mension, as shown in Table [T] These values agree with the dimensions given 
implicitly in |llj . 



2.3 Laplacian 

In |161 117j Laakso spaces are described as projective limits of quantum graphs 
and it is shown how to extend a compatible family of self-adjoint operators 
on the approximating quantum graphs to a self-adjoint operator on the limit 
space, i.e. the Laakso space. It was also shown how to use the spectrum with 
multiplicities of the operators on each quantum graph to determine the spectrum 
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3i 




2 


2 


3 


log(Q) 
log(i) 


[2,3,2,3,...] 


iog(24) 
log(6) 


[3,2,3,2,...] 


log(2i) 
log(6) 



Table 1: Hausdorff Dimension for Laakso Space associated with given sequence 
of : 



with multiplicities of the operator on the limit space. A quantum graph is a 
metric graph with a Hamiltonian operator, as described in [SI [TO], the simplest of 
which would be the Laplacian operator, i.e.: a Hamiltonian without a potential. 

On each metric graph, F n , consider the space of functions defined on the 
collection of edges each treated as a line segment. Define an operator on these 

,2 

function by A„ = — -gy. To make this a self-adjoint operator we need to also 
specify a suitable domain. A function is in Dom(A n ) if it is continuous every- 
where, continuously twice diffcrentiable on each line segment, and has Kirchhoff 
matching conditions at the nodes. Kirchhoff conditions require a function's with 
directional first derivatives summing to zero at nodes. 

Definition 2.2. A Laakso space, L, is a projective limit of the F n . Also there 
exist projection maps $„ : L — > F n for all n. Thus any function on F n can 
bepulled back to a function on L by writing f o $„ = / : L — > R. The pulling 
back is under the projections <£>„. 

By Theorem 7.1 in [17] those functions in Dom(A) that are pull backs are 
dense, so a complete set of eigenfunctions can be taken from this set. A con- 
sequence of this is that we can numerically approximate the spectrum of the 
Laplacian on the Laakso spaces working on some F n . Computations of these 



approximations are described in Section 3.2 along with calculations described 



in Section 3T Tables [2] and [3] show calculated values of the spectrum of the 



Laplacian for Laakso space. 



3 The Spectrum of A 

Theorem |3.1| gives the spectrum and associated multiplicities of the Laplacian 
operator by considering A„ on F n and on any Laakso space. We devote the rest 
of the section to proving the theorem. Following the analytic arguments are 
details of computational experiments carried out before the analytic results were 
available. We use an iterative, computer-assisted process to find the bottom end 
of the spectrum on a number of specific Laakso spaces. In all cases the computed 
results and analytic results agree within the precision of the computations. 
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Table 2: Calculated Values of the first 20 Eigenvalues for {j n }^Li = 
{2, 3, 2, 3, ...} with multiplicity, to, the iteration value, n, and the expected value, 
A. As n increases, the observed eigenvalues converge to the expected result. 



Theorem 3.1. Given any Laakso space, L, with associated sequence {ji}f =1 , 
the spectrum of A on Dom(A) is 

00 OO OO OO OO 

L>V}u u U{(*+i/2)V#u u U{ fcV/ «} 



k=0 n=l k=Q n=l k=l 

00 00 00 00 



u(JU{^} u UU 

n=2 k=l n=2 k=l 



(3.1) 



with associated multiplicities: 

1, 2", 2"- 1 (j n -2)/„_ 1 , 2"- 1 (/„_ 1 -l), 2"- 2 (/„_ 1 -l) (3.2) 
respectively. 

This does correct a typographical error in the similar statement given in |16j . 

3.1 Counts of Eigenvalues and Multiplicities 

In order to determine the spectrum of the Laplacian on the Laakso space, the 
approximating quantum graph is considered as a collection of simpler parts. In 
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Table 3: Calculated Values of the first 20 eigenvalues for given sequences of jVs 
with multiplicity, m and the expected value, A. 

A={2,3,2,3,...} B={3,2,3,2,...} C={3,4,3,4,...} D={4,3,4,3,...} 

[IH] it was determined that three distinct shapes with appropriate boundary 
conditions could be used to construct any quantum graph representations of 
Laakso spaces, save Fq, which is treated as a special case. Definition |3 . 1 1 defines 
these three shapes shown in Figure [5] with their respective boundary conditions 
which are forced by the Kirchhoff matching conditions and the orthogonality 
requirements that assign an eigenfunction to a given representation level. These 
orthogonality conditions were discussed in detail in [16] . In short they allow the 
counting arguments to count an eigenfunction only once. 

Definition 3.1. (a) A shape is a connected quantum sub-graph, as shown in 
Figure [5| In that figure the "D" denotes a Dirichlet boundary condition at that 
node and an "N" the Neumann condition. 

(b) A V is the shape consisting of three nodes: two nodes in a column and 
the third node a second. The two nodes in the first column are degree one and 
the node in the second column is degree two. When a V is in F n , it shares its 
degree two node with another shape thus making it a degree four node, as seen 
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V" "Loop" "Cross 



Figure 5: Constructions of v's, loops, and crosses along with associated bound- 
ary conditions 

in Figure^ 

(c) A loop is the shape that consists of two nodes each of degree two. The 
nodes are connected to each other by two cells, creating a loop. When a loop is 
in F n , both degree two nodes are shared as degree two nodes for another shape 
thus making them degree-four nodes, Figure^ 

(d) A cross is the shape consisting of six nodes four of degree two and two of 
degree four. The degree two nodes each have a cell connecting the node to each 
of the degree four nodes . Notice the cross is the only shape containing nodes 
of degree four in the subgraph. When a cross is in F n the degree two nodes are 
shared with degree two nodes of another shape thus making them degree four 
nodes, as in Figure^ 

Before determining the spectrum of A„ on the three shapes, we must first 
establish the following proposition which describes how these three shapes are 
involved in the construction of F n and L. 

Proposition 3.1. (a) Any node in any quantum graph approximating a Laakso 
space is either of degree one or degree four. 

(b) For any degree one node in F n , a V is produced in F n+ \. 

(c) For any degree four node in F n , a cross is produced in the construction of 
F n +i- 

(d) Any cell in F n produces j n +x — 2 loops in F n+ i between the V 's or crosses 
produced by the nodes in F n . 

(e) For n > the number of nodes in F n is N n — 2"~ 1 (/„ + 3). 

Proof. (a) A degree one node in gives rise to two degree one nodes in 

F n as an immediate consequence of the construction. Similarly a degree 
four node gives rise to a single degree four node in F n . The new nodes 
in F n that are not nodes in a copy of F n _i are the identification of two 
degree two nodes, hence of degree four. 

(b) In the construction of F n+ i, the cell connected to a degree one node is 
split into jn+i intervals by adding (j n +i — 1) nodes. Then the graph is 
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duplicated yielding two rows of cells connected between j„+i columns with 
two nodes in each, all of degree two except for the nodes at the end of the 
cell. The original and duplicated cells are connected at the newly added 
nodes. Thus the original node of degree one from F n remains degree one 
in F n+ i. The graph around the original node and it's duplicate is a "V." 

(c) In the construction of F„ +1 , the four cells connected to the degree four 
node in F n will be split into j n +i intervals. To construct F n+ i, F n is 
duplicated, new nodes inserted, and connected at the newly added nodes. 
Thus the original degree-four node remains degree four and is duplicated, 
creating two degree-four nodes. The graph around the original node and 
it's duplicate is a "cross." 

(d) Parts b and c account for two of the j n +i intervals. The rest produce 
loops. So, there are j n +i — 2 loops in F n+ \ for every cell in F n . 

(e) We will induct on n. The unit interval, Fq, has two nodes. Suppose that 
F n - X has iV„_x = 2"- 2 (/„_ 1 +3) nodes. Then N n = 2 x N n - 1 +2 n -\j n - 
l)7 n _i, the nodes from the two copies of -F n _i plus the new nodes of which 
there are j n — 1 new nodes per cell in F„_i and there are 2™~ 1 /„_i cells 
in F n _i. This simplifies to the claimed formula. 

□ 

We now generalize the results from |16| in three lemmas that give the eigen- 
values and multiplicites (counts) for each of the three shapes. 

Lemma 3.1. For any n > 1, the number of V's in F n is 2". The eigenvalues 
for this shape at this level are: 

{[/„(fc + l/2)7r] 2 : fc = 0,l,...}. (3.3) 

Proof. We prove the count by induction. F\ is constructed out of Fq, which is 
a single cell connecting two degree one nodes. This implies by Proposition |3.1| 
that F\ will have 2 V's. Now assume that for some arbitrary n > 1, the number 



of V's in F n is 2™. From Definition 3.1 the V is the only shape that has a 



degree one node. Furthermore, it has two degree one nodes. From Proposition 
3.1 we know that each degree one node in F n produces a V in F n+ \. From [16] 
the shapes defined in Definition [3T] are all the possible shapes in the graphs, so 
there cannot be any degree one nodes from any other shape. So the number of 
V's in F n+ i is twice the number of V's in F n . So F n+ i has 2 n+1 V's. 

In order to get the spectrum of A„ restricted to a V we look at the functions 
in this domain that are orthogonal to the functions expressible on an interval. 
These functions have the property that the values on the top branch are the 
negative of the values on the lower. We therefore need only consider the top 
branch, as it fully determines the behavior on the bottom branch. This top 
branch is one interval, and has Neumann boundary conditions at one end and 
Dirichlet boundary conditions at the other. The length of the cell is I" 1 . So we 
are looking for eigenfunctions on intervals of length l~ x with zero derivative at 
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one end and zero value at the other. These come in the form cos(/„(fc + 1/2)-kx) 
where k = 0, 1, . . . and x £ [0, J^" 1 ]. The eigenvalues in (3.3) are now obtained 
in the usual way. □ 



Lemma 3.2. For any n > 1, the number of loops in F n is 

2"- 1 (j„- 2) (/„_!). 
TTie eigenvalues for this shape at this level are 

{[Inkir] 2 :k=l,2,...}. 



(3.4) 



(3.5) 



Proof. By Proposition |3.1| every cell in F n _i produces j„ — 2 loops in F n . In 
order to know how many loops are in F n , the number of cells in F n —\ are counted 
and multiplied by j n — 2. The number of cells in F n were already counted in 
Proposition 2.1 and shown to be 2 n (I n ). Substituting in n — 1 for n in this 



expression and multiplying by j„ — 2 gives (3.4). 

In order to get the spectrum of A„ restricted to a loop we look at the 
functions in this domain that are orthogonal to the functions expressible on an 
interval. Again, these functions have the property that the vales on the top 
branch are the negative of those on the lower. As was the case with the V 
above, the orthogonality condition imposed on the functions reduces the ques- 
tion to only considering the top interval of length I~ x with Dirichlct boundary 
conditions. The cigenfunctions that fit these conditions are s in(/ n fc7rx) with 
k = 1,2,... and x £ [0, 1~ ]. These result in the set defined in (3.5). 

□ 

Lemma 3.3. For any n>2, the number of crosses in F n is 

2"- 2 (/„_ 1 - 1). (3.6) 
There are two sets of eigenvalues for this shape at this level. They are 

1 



:(I n kn) 



k = 1,2, 



with multiplicity one and 



with multiplicity two. 



{[I n kw] 2 :k= 1,2,...} 



(3.7) 



(3.8) 



Proof. From Proposition 3.1 crosses in F n appear only where there were degree 
four nodes in F n ^\. Therefore, to find the number of crosses in F n , we will 
count the number of degree four nodes in F n -\. By Proposition 3.1 every node 



in a quantum graph approximating a Laakso space is either of degree one or 
degree four. Therefore, subtracting the number of degree one nodes from the 
total number of nodes will give the number of degree four nodes. From the 
same proposition, the total number of nodes is 2™~ 1 (/„ + 3). We have seen 
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already that in F n degree one nodes only occur in V's and that for every v 



there are two degree one nodes. From Lemma 3.1 that there are 2 n V's in F n . 
Therefore, there arc 2" +1 degree one nodes in F n and 2 n ~ 1 (I n + 3) — 2™ +1 = 
2"~ 1 ( / « + 3 - 2 2 ) = 2"- 1 (/„ - 1) degree four nodes. Substituting n - 1 for n 



in this last expression gives (3.6). We note that this lemma is stated only for 
n > 2 because F± never has a cross since there are only degree one nodes in Fq. 

To obtain the spectrum of A„ restricted to the cross, we must consider the 
functions in the domain of the Laplacian on the cross. We can think of the cross 
as two X-shapes, (such as Fx in Figure [T]) connected at their four outer nodes. 
The orthogonality conditions from [TB] force the function on the bottom X to 
equal the negative of the function on the top X. The value of the function on the 
top of the X determines the value of the function on the bottom. The width of 
the X shape is 2/~ 1 and will have Dirichlet boundary conditions at the degree 
two nodes. Any function can be decomposed as symmetric and anti-symmetric 
with respect to the upper and lower branches of the X. We consider the two 
cases in turn. 

In the symmetric case, the function is the same along the top and bottom 
branches of the X. Therefore we need only to look at the top branch as it fully 
determines the bottom branch. Here we are looking for eigenfunctions on an 
interval of length 2/" 1 and zero at the boundaries. These are sm(^I n kirx) with 
k = 1, 2, . . . and x € [0, 2/~ 1 ]. The associated eigenvalues to these functions are 



those given in (3.7) 



In the antisymmetric case, the function value horizontally along the bottom 
branch of the X is the negative of the value along the top branch. At the 
central node, where the two branches meet, these two values must equal, so 
they must be zero. We then effectively have the X broken up into two V's of 
length J" 1 but with Dirichlet boundary conditions at either end. Looking at 
one of these V's, we still have the value along the bottom branch equal to the 
negative of the value along the top, so we consider only the top branch. Here 
we look for functions of length J" 1 with Dirichlet boundary conditions at both 
ends. This has already been done in Lemma |3.2| for the loop shape. There 
we got sm.(I n kirx) with k = 1,2, ... and x £ [0,/" 1 ] as the eigenfunctions and 
{[/„£:7r] 2 : k = 1,2, . . .} as the spectrum. This spectrum has multiplicity two 
because there are two halves of in the cross. □ 

Now we must consider the graph F and the eigenvalues it contributes to 
the spectrum. This graph is just the unit interval, and has Neumann boundary 
conditions forced by the Kerchoff matching conditions. So we are looking for 
eigenfunctions on intervals with length one and zero derivative at either end. 
These come in the form cos(fc7rcc) where k = 0, 1, . . .. This results in the following 
spectrum with multiplicity one: 

{[kir] 2 :fc = 0,l,...} (3.9) 

Table [4] we summarizes the results of these lemmas. In order to obtain 
the full spectrum with multiplicities, these sets must be combined with the 
multiplicities over all n > 0. Hence, Theorem |3 . 1 1 holds . 
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Shape 


Count 


Spectrum 


n Value 


F 


1 


{[tor]* :k = 0,1,...} 


n = 


V 


2 n 


{[/n(fc+|)7rf : fc = 0,l,...j 


n > 1 


Loop 


2"- 1 (j„-2)(/ n _ 1 ) 


{[I n k7vY : fc= 1,2,...} 


n > 1 


Cross 


2"- 2 (J„_ 1 -l) 


{[i(i„fc7r)f :fc = l,2,...} 
{[J„kf:t = l,2,...} x2 


n > 2 



Table 4: Summary of Lemmas 3.1 through 3.3 



3.2 Numerical Computations of the Spectrum 

A MATLAB script described in [IB] calculated the spectrum of the Laplacian 
for constant j n for all n € N by producing the incidence matricies of the ap- 
proximating graphs. We modified this script to handle general Laakso spaces. 
As in the original script, the eigenvalues are calculated using the eigs function, 
which is based on ARPACK (see users guide [IT)]). 

Quantum graphs approximating the Laakso space associated with the se- 
quence ji=[2,3,2,3,...] is shown in Figure [4] and the first twenty eigenvalues of 
the Laplacian on this space are shown in Table [2] The first twenty eigenvalues 
of Laplacians on other Laakso spaces are shown in Table [3j These computations 
agree with the calculations found in Lemmas |3.1| through |3.3| 



4 Heat Kernel 

Given a Laplacian on a Laakso space, the trace of its heat kernel will be obtained 
following the same procedure used in [5] where the heat kernel's trace was found 
for the diamond fractal. From [T5] the heat kernel of the Laplacian is 

p(t,x,y) = ^2 ^k,i(y)^k, m {x)e~ tEk 

k,l,m 

where ipk.i and ipk.m are L 2 — normalized eigenfunctions of A. The heat kernel 
on Laakso spaces will be further studied in |18j where continuity and bounds 
will be proved. The trace of the heat kernel Z(t) is defined in [2] as 

Z(t) = { p(t, x, x)dx = J2 9ke~ Ekt , (4.1) 
J k 

where are the eigenvalues of the Laplacian on the fractal and gk are the 
respective multiplicities associated with those eigenvalues. Associated with the 
heat kernel is the spectral zeta function also defined in [5] from the heat kernel 
as 

1 r°° 

C(s,7) = ^/ t*- x Z{t)e-*dt (4.2) 
r (s) Jo 
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where T(s) = J °° t s 1 e t dt is the gamma function. We set 7 = throughout 
the rest of this paper. Substitute (4.1) into (|4.2|) to obtain 



1 r°° 

r ( s ) ^0 

= F , g " ■ (4.3) 

The next step in any specific example is to simplify the spectral zeta function by 
recognizing Riemann zeta functions, Cr( s ) = — ji an d identifying the other 

n=0 n 

terms as geometric series. 

Definition 4.1. Define r = lim (Z^) 1 /™ iw/ien t/iis Zzmii exists. In the case of 

n— >oo 

self-similar spaces r is the contraction ratio since I^ 1 is the diameter of each 
cell. There is for any value of r a sequence j n that will produce that value. 

Once all of the series are simplfied, the poles of £(s, 0) in the complex plane 
can be calculated. The poles of £(s) for the diamond fractals are given in [5] as 

dh 2iirm „ 

s m = ~r + ti—> m e z ' ( 4 - 4 ) 

a w d w in r 

where dh and d s are the Hausdorff and walk dimensions respectively. 

Since the spectral zeta function was expressed as an integral of Z(t), applying 
an inverse Mellin transform [4 allows the heat kernel to be expressed as 

1 ra+ioo 

Zn(t) = — / C D {s)T{s)r s ds. (4.5) 

By the Residue Theorem, Zrj(t) is the sum of the residues of (D(s)T(s)t~ s . The 
residue must also be calculated at s — (a pole for T(s)) and at s = 1/2 (a pole 
for the Cr(2s) term in £(s, 7)). 

It is known that (r(s) = (r(s) and T(s) = F(s) for all complex s. Thus, 
the residues from s m and s_ r „ are complex conjugates of one another; therefore 
their sums equal twice the real part of the residue of s m . The complex values 
of the trace of the heat kernel yield oscillatory behavior in the heat kernel. We 
shall observe what happens to the heat kernel as t — > 0. Therefore, only the 
leading term with the most negative real power of t as well as any constants are 
included. For example, a result of [2] shows that for the diamond fractals, the 
trace of the heat kernel is 

Z D (t) ~ Cn(0) + ^-^-l^ (ao + 2Re( ai t-^ d ^)) + .... (4.6) 

This shows that the dominating power of t in the leading term as t — > is 
—d s /2. This incidentally is the complex dimension introduced in [12]. We shall 
now perform the same calculation for general Laakso spaces. 
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5 The Trace of the Heat Kernel on Laakso Spaces 

In Laakso spaces the dimensions and value of r are not always known by other 
means. Therefore the poles will be calculated analytically and the results will 
provide information about the Hausdorff, walk, and spectral dimensions. We 
now give the leading powers of the trace of the heat kernel for the Laplacian on 
a general Laakso space. 

Theorem 5.1. For the Laakso space associated to the sequence {ji}°Zi the trace 
of the heat kernel is 

oo oo 

Z(t) = ^2- 1 (/„_ 1 -l)^ e - fc2 ^ t 



OO 



k=l 



oo oo 



2"- 2 (/„-i - 1) ]T e~i k * W + ^ 2" ]T , 



-(fe+l/2) 2 7T 2 /„t 



n=2 
oo 



k=l 



+ 2»- 1 /„_ 1 ( 3 n ~ 2) £ e-**** + Y, < 
with an associated spectral zeta function 



n=i k=o 

oo 

_ — K 7T t 



(5.1) 



Us) 



00 2"- 1 (/ n _ 1 )(2 2s - 1 +j n - 1) + 2»- 1 (|2 2s - 3) 



E 



\n=2 

2 2 S + 1 _ 4 + • 

j ±L 

Jl 



J2s 



1 



(5.2) 



Proof. The spectrum of the Laplacian on various Laakso spaces was given in 
Table Has 

00 00 ( £,2^2 j2 ~| oo oo 



oo oo oo oo ^ 

-(Al)= UUW u UU - 

n=2fc=l n=2fc=l 



UU |J{(*+l/2)ViS} 

n=l fe=0 



oo oo 



uUU^ u U A2 



(5.3) 



n=l fc=l 



fc=0 



with respective multiplicities 

2"- 1 (/„_ 1 -l), 2"- 2 (/„_ 1 - 1), 2™, 2"- 1 / n _ 1 (j„-2), 1. (5.4) 



Direct substitution of the these values into (|4.1|) gives the heat kernel. By 
the associated sp© 

oo oo 

&(*)= EE 



(4.3) the associated spectral zeta function is 

2 n - 1 (/„-i - 1) + 2"- 2+2s (/„_ 1 - 1) 



n=2 k=l 

oo oo 

+EE 



(7 2 fc 2 7T 2 ) s 



2^+25 



(n(2k + l) 2 ir 2 ) s ' (I 2 k 2 ir 2 Y 



oo oo 



EE 



2"- 1 /„-i(j„-2) 



oo 1 

V 1 . 

k~l ( fc27r2 ) s 



(5.5) 
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This can be simplified by identifying Riemann zeta functions and, in certain 
cases, a Dirichlet Lambda function [T]. Then the function can be manipulated 
in to one sum. 



Us) 



00 2»- 1 (/„_ 1 )(2 2s - 1 + j n - 1) + 2"-!(l2 2s - 3) 



E 



+ 



\n=2 

2 2s + l _ 4 + - x 

72s 
Jl 



J2s 



as claimed in the proposition. 



(5.6) 



□ 



This expression of the spectral zeta function cannot be simplified further for 
a general Laakso space with an arbitrary sequence of jVs. However, it provides a 
common starting point for Laakso spaces in which the sequence of jVs is known. 
The next step is to locate the poles of the spectral zeta function. Recall that 
only poles which yield the most negative real power of t are considered since 
they produce the dominating behavior as t — > in the trace of the heat kernel. 
There are in general more poles than these. 

Proposition 5.1. Of all the poles of the spectral zeta function in the complex 
plane, the poles that yield the most negative real power of t in the leading term 
of the trace of the heat kernel are located at 



log 2r + 2?:im 
log r 2 



(5.7) 



for integer m. 



Proof. The trace of the heat kernel requires the most negative power of t which 
corresponds to the poles with the greatest real component due to the t~ s term 



in (4.5). The proof of the proposition relies on analyzing the series in (5.2) 



to find the poles. Note that the in the series in (5.2), the numerator has two 



terms; values of s will be calculated that will make the denominator grow at 
the same rate as the numerator. The value of s with the greatest real part are 



the poles that will be used. First rewrite the denominator of (5.2) as I^-I^ 
where 2s = s\ + S2- Select s to match the rate of growth for the first term in 
the numerator. To make I^ 1 grow at the same rate as I n —i, s i should equal 1. 
To have J* 2 grow at the same rate as 2" _1 S2 should equal log r 2. Therefore, for 
the first term, 2s = 1 + log r 2. It can be verified that any poles from the second 



term in (5.2) will not have a real component as large as this pole. Therefore, 



the real part of the poles that will yield the desired leading term in the trace 

log2 — 2 2 . Including 2mm in the numerator 



of the heat kernel are s = | 



2 log r 



gives all of the complex values of this pole 



□ 



Corollary 5.1. The dominating t term in the trace of the heat kernel has power 
—ds/2 = ~Re{s m ). 
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Proof. Recall from Section [4] that the dominating power of t in the trace of the 
heat kernel as t — > was —d s /2. A result of Proposition 5.1 is that the greatest 
real component of the poles yield the dominating power of t. Therefore, when 
calculating the residue of Ci( s )r(s)i~ s at the pole obtained from (5.7), the real 
power of t will be precisely —Re(s m ). But as stated at the beginning of the 
proof, it is also equal to —d s /2. Thus, d s /2 — Re(s m ) □ 

Corollary 5.2. The spectral dimension of any Laakso space with {ji} such that 
r exists is 

d s = (5.8) 
logr 

The walk dimension, d w , is 2 for any Laakso space which implies dh = d s . 



Proof. This follows directly from Proposition 5.1 and Corollary 5.1 since d s /2 

log 2r 



Re(s m ) = j§||£ which implies d 



og 2t r-Q^ wa vjj dimension is a result in 

log r 



It does agree with d% and d s via the Einstein relation 2dh/d w = d s 



□ 



The next two subsections give the trace of the heat kernel for two specific 
Laakso spaces where the sum in the spectral zeta function can be evaluated 
exactly: j = 2, and j = {2, 3, 2, 3...}. 

5.1 Laakso Space with j=2 

Proposition 5.2. For the Laakso space where at each iteration j — 2, written 
L2, the trace of the heat kernel is 



Cl 2 (0) 



1 + 2Re 



6Cfl(2+^)r(i+ I |p [ ) 



16t log 2 

Proof. For Laakso spaces with a fixed j = 2, /„ 
(5.1) and (5.2) gives the following two equations 

00 00 

Z L, (*) - E [ 22 "~ 2 ~ 2 " _1 ] E e_t22 " fe; 
n=l fc=l 
00 00 

+ [2 2n ~ 3 - 2"- 2 ] 

n=2 

00 00 



2™. Substituting these into 



k=l 



n=l 



k=0 



and 

Cl 2 M) 



( R (2s) /4(2 2s -! + 1) , 6(2 



r 2s 



is2-l 



i2s 



4 s (4 2 -4) 4 s (4 s -2) 



which has poles 



27rim\ 
log 4/ 



27rim\ 
log 4 J 



, Vm S Z. 



(5.9) 



(5.10) 



(5.11) 
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6 



Figure 6: Heat kernel Zl 2 , normalized by the leading non-oscillating term for 
the j — 2 Laakso space. The variable s is on the horizontal axis. 



By Proposition |5.1[ only the second term of the above equation with the 
greatest real part contributes to the leading t term. Then an inverse Mellin 
transform is applied just as in Theorem |5.1| Table [5] shows the residues of the 



integrand after the inverse Mellin transform at the poles given in (5.11 ) as well 
as s = from the T(s) term. Again we add complex conjugates and take only 
the most negative powers of t. Notice that when adding the residue from the 



poles in (5.11) only the poles with real part one contribute to the heat kernel, 
the others with the exception of zero and one half have residue zero. Once all of 
the residues are simplified we obtain the expression for the trace heat kernel's 
leading term as given in the proposition. □ 

The power of 1 of t in the denominator of the leading term in the proposition 
implies that the spectral dimension d s for L2 is 2. Knowing that the Hausdorff 
dimension dh = 2 from Tabic [lj we conclude that the walk dimension d w — 
2 since d s = 2dh/d w . Notice the poles and dimensions of this fractal were 
explicitly calculated. But Corollary |5.2| and Proposition |5.1| yield the same 
result once r = lim J, 1 /™ = lim (2 n ) /n = 2 is known. 
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Residue 







1/2 


8\/7fVt 


1 


16tlog2 


i _i_ 2irim 


i 6C H (2+^)r(i+^) 


log4 




1 _j_ 'ZlTT 

2 ~* log 4 


for m 7^ 



Table 5: Residues of the integrand of the inverse Mellin transfrom for given 
poles of the spectral zeta function for Laakso spaces with a fixed j = 2 



5.2 j={2,3,2,3...} Laakso Space 

Proposition 5.3. For the Laakso space with jik — 3 and j2fc-i = 2 where 
k > 1, the trace of the heat kernel is 

Z L (t) ~<l(0)+ (5.12) 

/p f 1 , log(2) \ f (a , 2 log(2) \ 
1 / 1 ^2 + log(6)J + iog(6) J 



1 I l"g(2) 

24t2 + io g(6) log ( 6 ) 



i+- 



r 



/ 1 , log(2) , 2irim \ a /, . 2 log(2) . iwim \ ' 
^2 log(6) "+" log(6) y V lo s( 6 ) lo S< 6 ) / 

, , 2 1og(2) 4„m 2-n-im 

77 l°g(6) + log(6)£ 155(61 



r, I 4 log(2) , Wm 1 | 2 1og(2) 4Tim > 

log(6) + log(6) _|_ 10- 2 lo S< 6 > iogf 6 ) + 12 

, , 2 1og(2) 4lrim 
2 i+ log(6) + log(6) 



Proof. In this case r = lim /V" = H m ( 2 ™/ 2 3"/ 2 ) 1 /" = \/6. This locates the 



n— foo 71— >oo 



poles with largest real part at s m = \ + '"^^ + ^j^p ■ 

The next step is to obtain and simplify the spectral zeta function associated 



with the trace of heat kernel given in (5.2). Since ji alternates between 2 and 



3, the following values can be directly substituted in, for any k we have: 
h=2 7 2 = 6 I 2k -2 = 6 k - 1 / 2 n = 2x6" I 2k = 6 k 



In preparation for substituting these values into (5.2 1 the sum is split into 
two sums, one over even n and the other over odd 
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(/„_! - 1) + 2"- 2+2s (/ n _! - 1) + T 



2" + 2"- 1 / n _i(j, l -2) 



_2(/ 1 -l) + 2 2s (Ji-l) + 2 2 



i 2s 



2 2 + 2h 



[2s 
1 2 



fc=2 



2 (/ 2fe -2-l) + 2 



2fc-3+2s 



-1) 



-,2fc— l+2s r,2fc — 1 



r2s 

d 2k-l 



+ 



2* k {I 2 k - 1) + 2 2fc - 2+2s (/ 2fc -i - 1) + 2 



2fc+2s 



_ 2 2fc +2 2fe-l 



(Jn-l) 



r2s 
1 2k 



(5.13) 



Substituting the known values of l2k, hk+i, 32k, an d J2k+i we obtain 



aw 



-2s 



■>2s 



„2s-l 



1 



12 

3 • 2 2s - 6 



24 



2 2s J 6 2s - 24 
6-2 2s -12\ I 



2 2.s 



(5.14) 



To apply the inverse Mellin transform as was done in (4.5| to obtain an 



expression for (t) . Then (t) can be calculated using the sum of the residues 



of Ci(s)r(s)i _s using the poles obtained in (5.7) as well as s — 1/2 (obtained 
from the Riemann zeta function, (r(2s) and s = (from T(s)). Table [6] lists 
the residues for those poles. 



1 , log(2) . 2-Kim 



log(6) 

'2 ~ I 



log(6) 1 log(6) 



Residue 



a(o) 



2 



^^2^ + 10 • 2*™ + 12) 

r(s fc )Cn(2s fc ) 3 4^-4 
t a fc7r 23 fc 81og(6) 4 3 fc 



Table 6: Residues of the integrand of the inverse Mellin transfrom for given 
poles of the spectral zeta function for the {2,3,2,3...}Laakso spaces 

The sum of the residues in Table [6] expresses the trace of the heat kernel 
with the leading terms as shown in the statement of the proposition. Note that 
the general terms in Zl (t) are shown in the proposition to indicate the behavior 
of the trace of the heat kernel in more detail. □ 

Corollary 5.3. The exponent oft, ^f^g , in the leading term as t goes to- 
wards zero implies that the spectral dimension of this Laakso space is d s — 
log 24/ log 6. The Hausdorff dimension for this fractal is given in Table [I] as 
d h = log 24/ log 6. 



Again, the same results are obtained by simply applying Corollary |5.2| and 
Proposition |5.1[ 



21 



Acknowledgments 



The authors would like to thank Alexander Teplyaev, Luke Rogers, Robert 
Strichartz, Shotaro Makisumi, Grace Stadnyk, Jun Kigami, and Naotaka Ka- 
jino. 

References 

[1] M. Abramowitz and I. Stegun, Handbook of Mathematical Functions: with For- 
mulas, Graphs, and Mathematical Tables, Dover Publications, 1965. 

[2] E. Akkermans, A. Comtet, J. Desbois, G. Montambaux, and C Texier, Spectral 
Determinant on Quantum Graphs, Annals of Physics, 284 (2000) 10-51. 

[3] A. Allan, M. Barany, and R. S. Strichartz, Spectral operators on the Sierpinski 
gasket I, Complex Variables and Elliptic Equations 54 (2009) 521-543. 

[4] G. Arfken and H. Weber, Mathematical Methods for Physicists, Elsevier Academic 
Press, 2005 

[5] M. Barlow and S. Evans. Markov processes on vermiculated spaces, in: Random 
walks and geometry, (ed. V. Kaimanovich), de Gruyter, Berlin, 2004. 

[6] D. Ford and B. Steinhurst, Vibration Spectra of the m-Tree Fractal, To appear 
in Fractals (2010). 

[7] K. Igudesman, Lacunary Self-similar Fractal Sets and Intersection of Cantor Sets, 
Lobachevskii Journal of Mathematics 12 (2003) 41-50. 

[8] N. Kajino, Spectral asymptotics for Laplacians on self-similar sets, J. of Func. 
Anal. 258 (2010) 1310-1360. 

[9] P. Kuchment, Quantum graphs I. Some basic structures, Waves in random media 
14 (2004) S107-S128. 

[10] P. Kuchment, Quantum graphs II. Some spectral properties of quantum and 
combinatorial graphs. J. Phys. A. 38 (2005) 4887-4900. 

[11] T.J. Laakso. Ahlfors Q-Regular Spaces With Arbitrary Q > 1 Admitting Weak 
Poincarc Inequality, Geom. Funct. Anal. 10 (2000) 111-123. 

[12] M. Lapidus and M. van Frankenhuysen, Fractal Geometry and Number Theory: 
Complex Dimensions of Fractal Strings and Zeros of Zeta Functions, Birckhauser, 
Boston, 2000. 

[13] M. Lapidus and M. Pang, Eigenfunctions of the Koch Snowflate Drum, Comm. 
in Math. Phys. 172 (1995) 359-376. 

[14] M. Lapidus, J.W. Neuberger, R.J. Renka, Snowflake harmonics and computer 
graphics: numerical computation of spectra on fractal drums, Internat. J. Bifur. 
Chaos Appl. Sci. Engrg. 6 (1996) 1185-1210. 



22 



[15] Lehoucq, R.B., D.C. Sorensen, and C. Yang, ARPACK Users' Guide: Solution 
of Large-Scale Eigenvalue Problems with Implicitly Restarted Arnoldi Methods, 
SIAM Publications, Philadelphia, 1998. 



[16] K. Romeo, B. Steinhurst, Eigenmodes of a Laplacian on some Laakso Spaces, 
Complex Variables and Elliptic Equations 54 (2009) 623-637. 

[17] B. Steinhurst, Dirichlet Forms on Laakso and Barlow Evans Fractals of Arbitrary 
Dimension arXiv:0811:1378. 

[18] B. Steinhurst, Uniqueness of Brownian Motion on Laakso Spaces, in preparation. 

[19] R. Strichartz, Differential Equations on Fractals: A Tutorial, Princeton Univer- 
sity Press, Princeton, New Jersey 2005. 



23 



